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3 ■ ABSTRACT 

in 



We explore the potential cumulative energy production of stellar mass black 
holes in early galaxies. Stellar mass black holes may accrete substantially from 
^ I the higher density interstellar media of primordial galaxies, and their energy re- 

lease would be distributed more uniformly over the galaxy, perhaps providing a 
different mode of energy feedback into young galaxies than central supermassive 
black holes. We construct a model for the production and growth of stellar mass 



O 

^ ■ black holes over the first few gigayears of a young galaxy. With the simplifying 

j^ ■ assumption of a constant density of the ISM, n ~ 10^ — 10^ cm~^, we estimate 

the number of accreting stellar mass black holes to be ~ 10^ and the poten- 
tial energy production to be as high as lO^^ergs over several billion years. For 
ly-^ ■ densities less than 10^ cm~^, stellar mass black holes are unlikely to reach their 

^ ' Eddington limit luminosities. The framework we present could be incorporated 

cn ' in numerical simulations to compute the feedback from stellar-mass black holes 

C^ ■ with inhomogeneous, evolving interstellar media. 

Keywords: Galaxies: evolution - Galaxies: stellar content - Stars: massive - 
Black hole physics - Accretion, accretion disks - Radiative transfer 



^ . 1. Introduction 

Most research on the early formation and growth of black holes has focused on the 
important issues of the growth and feedback of intermediate and supermassive black holes; 
relatively little attention has been paid to the effects of stellar-mass black holes and their 
potential energy production. Estimates of the number of stellar mass black holes presentl y 



in the Galaxy are as high as 10^ to 10^ (jAgol fc Kamionkowskil 120021 : IChisholm et al.ll2003l ) 



suggesting that the nearest is only parsecs away, but accreting feebly from the modern 
interstellar medium. The cumulative energy production of 10^ IOMq black holes could rival 
that of a single central IO^Mq black hole. 
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The first stars that formed were extremely low metallicity and arguably high mass 
( iBromm et al.ll2002l ) Population III stars. These stars may have only produced intermediate- 
mass or supermassive black holes, if any, but turbulence may broaden the initia l mass func- 



tion s o that lower mass stars and black holes form even in the first generation (j Clark et al. 



2011uPrieto et al.ll201ll ). The next generation of stars. Population II, were slightly higher in 
metallicity and would form stellar mas s black holes. This epoch might have begun at redshift 
in excess of 15 (JBromm fc Loebll2006l ) during the era of reionization, when the interstellar 
medium number density, Uism, niight still be large. In the early formation of these galaxies, 
their interstellar mediu m (ISM) density was probably much higher than it is currently, per- 
haps riism ~ lO^cm"^ (JGreif et al.ll2008l ). With ISM densities this high, stellar- mass black 
holes could potentially accrete substantial amounts of matter over the course of the early 
evolution of galaxies and emit energy back into the galaxy and intergalactic medium. Since 
the number of stellar mass black holes could be as high as 10®, their cumulative energy pro- 
duction could be extremely high, given a favorable accretion and feedback rate. Since stellar 
mass black holes would be distributed more evenly over the galaxy than a single supermas- 
sive black hole, the energy feedback may impact the evolution of a galaxy in different ways 
than would a singular source. Stellar- mass black holes also deliver feedback in a different 
mode than supernovae, the energy of which is primarily kinetic, and, unlike supernovae, 
black holes persist. The rate of star formatio r i in galaxies is about a factor of 100 less than 
provided by basic estimates (JKennicuttl Il998l : ISomerville et al.l l2008l ) . This suggests a feed- 
back process that has not yet been identified. The process we outline here may contribute 
to this regulation of star formation in early galaxy evolution. 



We a pply the prescrip t ions o f black hole accretion with feedback from (JMilosavljevic et al 



20091 ) and iPark fc Ricottil (1201 ll ) to the stellar mass black hole case to approximate the en- 
ergy production in early galaxies. With these estimations, we find that the potential energy 
production of the stellar mass case could be on the order of 10^" ergs over ~ 2 Gyr. Section 
2 gives an outline of our model for the birth and growth by accretion of stellar mass black 
holes. Section 3 addresses the rate of accretion of stellar mass black holes including the ef- 
fects of radiative feedback and gives estimates of the number of black holes, their luminosity, 
and their cumulative energy liberated as a function of time. Section 4 discusses the possible 
constraints of the Eddington limit and §5 gives our conclusions. 
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2. Star Formation Rate and the Black Hole Birth Rate 

The total number of stars born per unit time per unit galactic mass at time, t, is given 
by: 

/■oo 
nit)= MMMS,t)^MMs)dMMS: (1) 

Jo 
where ri^{MMs,t) is the rate of birth of stars per unit galactic mass at a given time with 
main sequence mass, Mms, and ^{Mms) is the mass function of main sequence stars with 
mass between Mms and Mms + dMMs- The observable star formation rate per unit galaxy 
mass, s^, at a given time is then given by: 

s,(t) = / ri,MMsHMMs)dMMs, (2) 

Jo 

where s^(t) has units of solar masses per year per unit galaxy mass for Mms measured in 
solar masses. 

Assigning a probability, Pbh{Mms), that a star of main sequence mass, Mms, forms 
a black hole, the black hole birth rate, rBH{t), and mass formation rate, SBH{t), per unit 
galaxy mass are given by: 

/■oo 
TBHit) = / ri,^iMMs)PBHiMMs)dMMS, (3) 

and 

/■oo 
SBHit) = / ri,^MMs)MBH,oiMMs)PBHiMMs)dMMS, (4) 

Jo 
where MBH,o{MMSiii Z) is the initial mass of the black hole that can, in principle, be a 
function of the main sequence mass of its progenitor, the epoch it was born, the metallicity 
of the progenitor star and perhaps other p arameters. Black holes in binary systern s have 



measured masses in the range 6-20 M© (jCasaresI l2006l : IShaposhnikov et al.l 120111 ) . For 



simplicity, we will scale our results with the assumption that all black holes are born with 
a single mass, Mbh,o = 10 Mq. We can then define the efficiency, bbh, to make black holes 
from a generation of stars at time t as: 

esH = 7-^ ^ Mbh,o 7-^, 5 



and thus 



rsHit) ^ -^eBuit). (6) 

The parameter eBH is proportional to the fraction of all stars that are 0-type stars, /o*, and 
the fraction of 0-type stars that eventually turn into black holes, Jbh- This efficiency could 
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vary with epoch and other parameters, but, again for simphcity and illustration we will take 
it to be constant. We adopt esH ~ fo^fsH ~ 10^'^. This corresponds approximately to 
making black holes from stars with Mms > 3OM0(that is, Pbh{Mms) being a step function 
at 30 Mq) for a Salpeter mass function with minimum mass of ~ 0.1 Mq. We thus find for 
the rate of birth of black holes per unit galaxy mass, 

rBH{t) ^ 1 X 10-S,(t)eBH,-3M^^_o,o yr'' (7) 

where Mbh,o.io is the mass of the initial black holes in units of 10 M0 s*(t) is in MQyr~^ per 
unit galaxy mass and eBH,-3 is the black hole production efficiency in units of 10^^. Below 
we will want to address the rate of star formation and the associated rate of production of 
black holes over an entire galaxy. We define RBuif) = MgairBuif) and S^{t) = MgaiSi,(t) so 
that 

RBH{t) ^ 1 X 10-^S,{t)eBH,-3M^l,,^,, yr-^ (8) 

where S^{t) is now the star formation rate integrated over a whole galaxy, expressed in units 
of solar masses per year. 



3. Accretion onto Stellar Mass Black Holes in Early Galaxies 



After stellar-mass black holes are born, they may grow by accretion from the early, dense, 
ISM. This accretion will be subject to feedback eff ects (JAlvarez et al.ll2009l : iMilosavljevic et al, 



2OO9I : lOstriker et al.l l2010l : iPark fc Ricottil I2OIII ) and the accretion rate will change as the 



black holes grow in mass. As mentioned in the Introduction, the density of the ambient ISM, 
rii^m,, could have been much higher in early galaxies, with some estimates as high as lO^cm 



flGreif et al.ll2008l : lMilosavlievic et al.ll2009l : IPark fc Ricottil l20 1 ll : IPebuhr et al.ll20ld . l2011f ). 
The ambient density is, in turn, a principal factor controlling the accretion. 

During the lifetime of a massive main sequence star, stellar winds will evacuate the local 
interstellar medium in the vicinity of the star. The progenitor star may explode rather than 
collapsing quietly to form a black hole. Once the progenitor reaches the end of its lifetime and 
forms a black hole, there will be some time, t^eiay, as the ISM refills the excavated volume, 
the bubble "pops" or the black hole drifts from its natal site. After this delay time, the black 
hole begins accreting from the ambient ISM. The time span of the delay may be diminished 
by neighboring stars exploding as supernovae and displacing the ambient gas and dust near 
the vicinity of the blac k hole. Estimates give t^eiay ~ lO^yr (JGreif et al.ll2008l : IWise fc Abel 
2008uGreif et al.ll2010l ). The physics governing this delay time is complex, but for simplicity 
we will take the delay time to be constant. Over the time spans we will consider, gigayears, 
this may not be an important effect, but we will formally keep the term through much of 
the analysis. 
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Recent computational models have presented accretion rates and the associated en- 
ergy production with radiative energy and momentum feedback for intermediate and su- 
permassive black ho l es accreting in ISM with the high densities exp e cted in young galaxie s 



( lAlvarez et al 



2009[ iMilosayiieyic et all l2009l : lOstriker et all l2010l : iPark fc Ricottil 120111). 
The models of iPark fc Ricottil (120111 ) are especially useful for our current purpose by pro- 
viding a convenient pa rameterization of the effects of f eedback as a ra t io wi th respect to 
standard Bondi/Hoyle (JBondi fc Hoylelll944| ) accretion. iPark &: Ricottil (l201ll ) give the ac- 
cretion rate for a stationary black hole of mass Mbh accreting from an isothermal gas of 
temperature Tism and density pism as: 



Mbh,o =< Kad >o Mb. 



isofli 



(9) 



where the radiative feedback efficiency factor, < Xrad >o, is the mean ratio of the actual 
accretion rate to the Bondi/Hoyle accretion rate for the same ambient temperature and 
density, and where the subscript refers to the assumption that t he black hole has velocity 
V = with respect to the ambient medium. iPark fc Ricottil (120111 ) normalize their results to 
the Bondi/Hoyle accretion rate for a stationary black hole in an isothermal gas with sound 
speed Cg that they take to be: 



M, 



BAso,0 



XbG MBijmpnism 



(10) 



wher e A^ is a coefficient of order unity. Taking A^ = e^/^vr 
20091 ). Equation [To] can be written in convenient units as. 



14.08 ( JMilosavljevic et al, 



MB,iso,o = 1-4 X 10-^M2on,,^,5T,,^/^MQyr-^ 



:ii) 



where Miq is the black hole mass in units of IOMq, nism,5 is the ambient particle density 
in units of 10^ cm"'^, and Tism,4, is the ambient temperature in units of lO^K, were we have 
taken Tism = mpcl/k. For a black hole moving with a velocity, v, one would have for the 
Bondi/Hoyle accretion rate: 



M, 



B,iso 



M, 



B,iso,0 



(t;2 + c2)3/2- 



(12) 



The results of iPark &: Ricottil (120111 ) can be scaled to results appropriate for a moving black 
hole by an appropriate scaling of the radiative feedback efficiency factor and the Bondi/Hoyle 
accretion rate as indicated below. 



Remarkably, iPark fc Ricottil ( l201ll ) (see also iMilosavljevic et al.l |2009|) find that the 
radiative feedback efficiency parameter, < X^ad >o, is independent of the the accretion 
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efficiency, r], essentially independenent of the mass of the black hole and insensitive to the 
ambient density. They give: 

< A,,, >o~ 0.047;;^/^ ( s x^W ) ■ ""^^^ ^ lO^cm-3, (13) 

and 

< A.,, >o~ OMTr'J', ( s x^IO^k ) ""^^"'^ '■ ''''"' ^ 10^cm-^ (14) 

where Tin is the time-averaged temperature at the accretion radius. There is an implicit 
weak dependence of < Xrad >o on the mass of the black hole through the parameter, Tj„. 
The average temperature in the inner HII region will be a function of the spectral index of 
the radiation with harder spectral indices giving s maller accretion radii and higher Tj„. We 



are considering smaller mass black holes than did iPark fc Ricottil (120111 ) for which the disk 



radiation will be harder. To a l low fo r this, we take T^ corresponding to the hardest spectra 



considered by lPark fc Ricottil (1201 ll ). with spectral index a = 0.5 in their Figure 9 and thus 
Tin ~ lO^K. We ignore any other scaling of < Xrad >o with Mbh, although such dependence 
may exist. With this scaling, we have: 

< A,,, >o~ 0.016T-^';^7;;J, (15) 

for riism > lO^cm"^, and a multiphcative factor of n^/^g for Uism < lO^cm"^. 



In the final version of their paper, jPark fc Ricottil (120111 ) revised Equations [T3I and 



Ml reducing the overall efficiency from 0.04 to 0.03 and the characteristic value of Tj„ to 
~ 4 X 10^ K. These changes would alter the results here by making the overall accretion 
efficiency greater by a factor of ~ 30. Because the value of Tj„ for the case we are considering 
is uncertain, we have preserved our normalization to T^ = 10^ K. Our results can be scaled 
to other values of Tj„ in the manner we present, but our result for the luminosity from the 
ensemble of stellar mass black holes may be underestimated. The energy is somewhat less 
sensitive than the luminosity. The time, too, for a black hole to grow to infinite mass might 
be smaller and more comparable to the initial accretion delay time, t^, that we have generally 
neglected. 

For the case of a moving black hole, the velocity correction to the accreti on rate as given 



i n Equ ation Il2l means that for given Tism and Cg, the accretion rate given by lPark fc Ricotti 



( 120111 ) is too large. The actual accretion rate for a moving black hole would be less, as if the 
accretion were occuring at higher effective values of Tism and Cg. For a given Tj^^, a moving 
black hole would accrete at a rate corresponding to a higher effective ambient temperature 
of a stationary black hole such that: 

t)2 _)_ ^2 

-'■ism,v=0 -'■ism,v>0 7^ ; l-'-^/ 

C 
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where Tism,v>o is the actual ambient temperature and Tisrn,v=o is the ambient temperature 
(a nd corresponding soun d speed) that would give the proper accretion rate in the calculations 
of iPark fc Ricottil ( 120111 ) for stationary black holes. If, for a given ambient temperature, the 
accretion rate for a moving black hole corresponds to accretion at effectively a higher ambient 
temperatur e for a stationary blac k hole, then the values of the radiative feedback parameter 
derived by iPark fc Ricottil (120111 ) for a given ambient temperature must also be scaled to 
that same, higher effective temperature or, from Equation [151 



< \-ad{Tism,V > 0) > — < Xrad >C 



v'^ + cl 



-5/2 



(17) 



where < Xj-ad >o is the radiative fee dback efficiency param eter for the same value of Tj^^, 
but for a stationary black hole from iPark fc Ricottil (l201l[ ). We can thus write for a black 
hole moving with a velocity, v: 



MbH =< Xrad{Tism,V > 0) > MB,iso =< Xrad >0 Mb. 



iso,0 



2 1 2\ -4 



(1^ 



The results of (jPark &: Ricottil |20 111 ) for a given ambient temperature are thus equivalent to 
those of a somewhat higher "effective" ambient temperature for a moving black hole. This 
correction factor is quantitatively important, but not qualitatively important for v ^ Cg. The 
steep dependence on this factor may become significant if the black holes move supersonically 
through the ambient medium. 

With Equations [TOl HH [13] and [HI we can write: 



M, 



BH 



FM 



BH-i 



(19) 



where the factor, F, is: 



F —< Xrad >o 



v'^ + c 



^^ e^/^7rG^m„n. 



p"'ism 



(20) 



or 



F - 2 2 X lO^^^ff^^s^^T"^ T^ n^/^ 



3/2 



2 1 2N -4 
V^ + cf 



1 -llisinjS 



<1, 



(21) 



with Cg = kTisml'fnp and the factor of n^^^^ would be nj^^.s for riismfi > 1- The coefficient, 
F, is a function of the ambient medium, but independent of Mbh (except for the implicit 
dependence through Tjn^s). The solution to Equation [19] is thus: 



MBHit) 



M 



BHfi 



[l-FMBH,o{t-tacc)Y 



(22) 



where Mbh,o is the initial mass of the black hole and tacc is the time when the black hole 
began to accrete. For a black hole born at time t^ for which accretion was delayed by a time 
tdeiayy then tacc = ^fe + t delay Figure 1 gives a schematic diagram of the growth of a black 
hole seed in this framework and Figure 2 gives the black hole mass as a function of time for 
various choices of the interstellar density, with other parameters held constant. 

The timescale for a seed black hole to grow to infinite mass is given (for nj^^, 5 < 1) by: 

1 / 2 I 2 \ ^ 

Note that this time scale increases very steeply with the velocity of the black hole. Writing 
the velocity correction factor as f^ = {v^ + cj.)/c1 and combining Equations |22] and [231 we 
can write: 

V too y I i 7.2/4Gj/r 

for fiducial temperature parameters, and for nism,5 < 1. In our model, the most massive 
black holes will be those that began growing from seeds at a time, tdeiay after the first stellar 
mass black holes began to form at to, or tacc = to + t delay With seed black holes of mass 
Mbh.o = IOMq, the most massive black hole will have a mass of about 14M0 2 Gy (at z 
3) after the first stellar mass black holes formed and about 17Mq at 3Gy. The mass would 
formally become infinite at 7.2 Gy for the chosen parameters, but conditions, especially riism, 
would have changed by then and the accretion might be limited by the Eddington limit. For 
the timescales of interest here and nism.,5 < 1, this extreme growth is not of interest. We 
return in §l]to discuss the possible effect of growth to accretion at the Eddington limit. 

We can now evaluate the number of black holes that will have formed by a certain epoch, 
their luminosity at that epoch, and the energy they will have emitted by that epoch. 

The number of black holes of a specific mass, MsHit) at a time, t, will depend on the 
rate at which the seed black holes were born and their accretion history. Equation [22] can 
be inverted to give the interval of time that a black hole of given mass has been accreting, 

_MBH{t)-MBHO _ ( Mbho\ 

Atacc - t - tacc - FMbH^oMbhH) - ^^ [} ' M^) ) ' ^^^^ 

Since a black hole of given seed mass was born at a time tj, = tacc — t delay = t — ^tacc — t delay, 
we can then write the time when a black hole of mass MBH{t) was born as: 

t, = t-t^l^l- JJ^^) - tdelay (26) 
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from which we can write: 



^^'^ ^ (27) 



Note the sign change here with respect to Equation [19] that resuhs from the difference in 
taking the derivative with respect to t, holding tb constant versus taking the derivative with 
respect to t}, holding the currrent epoch, t, constant. The rate at which black holes of mass 
MBH{t) at time t were born is the rate at which their seeds were born at time, tb, which can 
thus be expressed as: 

RBHitb) = RBHyt — too I 1 7T 777 1 " tdelay)- (28) 

V MBH\tjJ 

Since black holes of larger current mass were born earlier, black holes with mass between 
Mbh and Mbh + ^Mbh were born between tb and tb - dtb- At current epoch, t, the number 
of black holes with mass between Mbh and Mbh + AMbh is thus given by: 

^^^^dMBH = '^iM^^^dp^dMBH = dJl^^J^dMBH, (29) 

ciMbh dAtacc uMbh dtb dtb dMBH 

where we have taken '^^^^^^f^"-^ = ^^^^^. With ^^^^^ = RBnitb), ^ = -1, and 
Equation [271 we can then write: 

^^^^^dMBH - RBHitb) j^^-^,. (30) 

The number of black holes in a galaxy with a given total star formation rate and hence 
total black hole formation rate born between the beginning of the epoch when stellar-mass 
black holes formed, to, and the epoch under consideration, t, is thus: 



NBHit)= / RBHit')dt'= / ^,/ uL dMBH+ / RBHit')dt', (31) 

Jto JM^Hn tlVlBHKt) Jt-t,„,„„ 



ft rMgH.maxit) 

RBH{t')dt' = / 

where tb{t') is given by Equation [26] and 



MBH,maxii) - Z 77 7 7 \T1~ ' *^^^) 

-L — (,t — Iq ~ i^delay)/i^c 



''oo 



is the maximum mass to which a black hole could have grown at time t, one which began 
accreting as soon as it could, a time tdeiay after the first stellar-mass black holes were born 
at to- The second integral over Rbh on the right hand side accounts for the black holes that 
were born during the last interval, t delay, before the epoch, t, under consideration that have 
not yet begun to accrete. The integral of Rbh over the full time interval, to to t, in Eqn [31] 
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yields the desired result for the number of black holes (especially in the trivial case for which 
Rbh = constant), but the first integral on the right hand side provides the framework for 
evaluating the luminosity and energy production, as described below. There is no additive 
term over the last interval of tdeiay in the computation of luminosity and energy, since these 
late-born black holes are, by assumption, not yet accreting. 



In principle, to evaluate the integral in Equation [3T] one needs to take into account the 
inhomogeneity and temporal variations in the ambient quantities that determine the factor, 
F, especially Uism and the velocity term, /^, and the temporal variation in the formation rate 
of black holes that itself is determined by variations in the star formation rate (Equation 
E]). Wh en young galaxies f requently collide, the star formation rate is expected to be very 
bursty (JDebuhr et al.ll201l[ ). To get a qualitative feel for the expectations of our framework, 
we will take the ambient conditions and the star formation rate to be constant in time. A 
representative star formation rate for a fiducial galax y of mass Mgai = 4 x 10^^ M© from 
the merging galaxy simulations of iDebuhr et al.l ( 1201 ll ) is S** ~ 10 M© yr~^ and hence from 



Equation [U Rbh ~ 1 x 10" 



NBH{t) 



,m; 



eBH,-3^''^BH,0,10 \ WMqvt- 



yr 



-1 



We can then write: 



R 



BH 



MBH.maxi't) 



dM^ 



BH 



F 

Rbh (t 

1 X lOhayr^BH, 



Mbh,o 

to) 



MBHity 



,M, 



(. 



s. 



(33) 



where in the last step we have taken t to be measured in Gyr and have neglected to compared 
to timescales of Gyr. We have taken the more cumbersome means to compute this quantity 
from Eqn. [HUfor illustration, but note that the result is consistent with the simple integration 
over the rate of production of black holes. We thus estimate that over a period of 2 Gy, a 
galaxy like the Milky Way could have produced of order 2 x 10^ black holes that had grown 
from seeds of 10 M© to various masses. Figure 3 gives the number of black holes as a function 
of time for various choices of the initial black hole mass, Mbh,o, with other parameters held 
constant. 

We can now use the same framework to estimate the luminosity of the ensemble of 
accreting stella r -mass black holes at the epoch, t. We will adopt the parameterization of 



Park fc Ricottil (120 111 ) with the assumption that there is sufficient angular momentum to 



form a disk near the black hole, so that disk-like efficiencies for turning m ass accretion rates 



i nto r a diated energy are applica ble. For critiques of this assumption, see iRuffert fc Arnett 



(I1994J ): iBeskin &: Karpovl (120051 ). We write for the luminosity of a single accreting black 
hole: 

L{t)=r,MBH{t)c\ (34) 
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with rj ^ 0.1. If a disk does not form, or forms only sporadically, the radiation efficiency 
would be correspondingly less than the fiducial value we assume here. 

The luminosity per unit galaxy mass radiated by all black holes with mass between 
Mbh and Mbh + dMsn at epoch, t, is given by: 



dLtotit) ^^^ _dNBH{t)^,^ ^o^,,2 



dMsHit) dMsHit) 

Using Equation |29l this can be written as: 

dLtotit) 



dMBH = —f^dMBHVc'FM'sj,{t). (35) 



dMBHit) 



dMsH = RBH{tb)vc'dMBH. (36) 



The total luminosity from all the accreting black holes born since to can then be obtained 
by integrating over all the current masses at epoch, t, to obtain: 

rMgH.maxit) 

Ltot{t) = / RBHih)vc'dMBH. (37) 

Once again, we can approximate the complex variation of the rate of production of black 
holes with a constant to obtain: 

Ltotit) ~ RBHTlc\MBH,ma.{t) - Mbh,o) = RBHr]^MBHfl f , ^ ' ^^ ' ^'^-^-V \ . (38) 

ytoo + to + idelay ~ t J 

Neglecting to and t^eiay^ taking t « too, and using Nbh ~ Rsnt gives a simple estimate of 
the luminosity of L^t ~ N BHif)'!]^^^ Mbh, o/ioo- Neglecting to and tdeiay we have for the total 
luminosity of a galaxy of constant star formation rate: 

T (+\ 8 X lO^^ergs-^Gy, 3/2 4 ( ^* \ /on^ 



where we have taken the timescale too from Equation [23] to be 7.2n^gJ^ 5/^-^5^,0,10 Gyr for 
fiducial temperature parameters and for nism,5 < 1- For t ~ 2 Gyr, the ensemble of about 
10^ black holes will produce ~ 10^^ erg s~^, as much luminosity as a modest Seyfert galaxy 
or a single black hole of about IO^Mq accreting near the Eddington limit. For comparison, 
supernovae provide an input of about Lsn ~ 10^^ erg yr~^ in a galaxy like the Milky Way 
for which the star formation rate is about 1 Mq yr^^, so a star formation rate of 10 Mg 







yj^-i jjiight give Lsn ~ 3 x 10^^ erg s~^. The black hole input at about 2Gy is comparable 
to, and might even slightly exceed the input power from supernovae. Figure 4 gives the 
luminosity versus time for various values of the interstellar density and Figure 5 illustrates 
the sensitivity of the luminosity to the velocity parameter, /^. 



-12- 



The total energy liberated by the number of black holes, Nbh-, accreting from to + 
is 

Etot = / Ltot{t')dt' = / RBH{tb)vc'dMBH. (40) 



'to+tdelay -J to+tdelay "^ ^'I BH .0 

We will again adopt the approximation of a constant value of Rbh to write 



io+tdelai 

Invoking 

FM'^^it'), (42) 



Etot ~ RbhVC^ I {MBH,maAt') " MBHfi)dt' . (41) 

•^to+tdelo'' 

dMBH,max{t') 

dt' 
we can write, 

'^'BHfi 

RbhVC f'--'"''BH,oi'-'o-tdeiay) MBH,max- Mbh,0,.. /.o\ 

J^totyt) = ^ / 772 (i^y^BH,max, (,4dj 

^ JMbh,o ^^^BH,max 

which becomes 

Etotit) = RBHVc'MBH,otoo S^-ln (l - ^~Y^'''°' ) + ^^^5_1^^| . (44) 

For to + 1 delay « t « too, this rcduces to 

Etot ~ 2RbhVC^Mbh,o (t + ^^ j ~ 2Ltotit)too ~ 2NBHit)r]c^MBH,o, (45) 

nearly independent of too- Again taking too from Equation [231 to be 7 .2n^^J^ ^f^M]^]j ^ ^q Gyr 
for fiducial temperatures and neglecting to and tdeiay on the timescales of interest, we have: 

Etot{t) ~ 4.1 X lO^Vgs RBHn^^^/sf,' 

X [-In (l - O.UtGyrn/^^J-^MBH,o,io) + O.UtGyrn,iij;;^MBHfi,io} .(46) 

For Rbh ~ 3.2 x 10-i^s-^eBH,-3MB^o,io ( iomoVO ' ^"^ ^'^'''' 

Etotit) ~ 1.3 X 10^°ergs eBH,-3MB^,o,io ( ioM^Jy,-i ) ""^ 

X [-In (l - Q.lAtGyrn]il^J;^MBHfi,io) + Q.UtGyr'4il^J-^MBH,o,io} -(47) 

For t ~ 2 Gyr after the beginning of stellar-mass black hole formation, corresponding to a 
redshift of about 3, we have 

Etotit) - 1.3 X lO^^ergs eBH,-3MBH,o,io ( ioM^Jy,-i ) ^.^m'^sfv 

X [-In (l - 0.28n^/^,5/-^MsH,o,io) + ^^M^H^J^^Mbh^q^iq] ■ (48) 
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For a density of rtism = 10^ cm~^, an initial black hole mass of 10 M© and taking the velocity 
factor, /„ to be unity, we get an energy liberated of Etotif) ~ 8 x lO^^ergs by 2 Gyr. This 
is equivalent to 8 x 10^ supernovae liberating about 10^^ ergs apiece, but the energy would 
be entirely in radiant, not kinetic energy, in our model. In practice, if the black holes are 
accreting by means of an accretion disk, as we assume, then some of the energy emitted may 
be in the form of jets and hence of kinetic energy. Figure 6 gives the total liberated energy 
versus time for various values of the interstellar density and other parameters held constant. 
Figure 7 gives the total energy for various choices of the velocity parameter, /^ at a density 
of nism,b = 1- Figure 8 gives the total energy as a function of the parameter, too- 



4. Eddington Limit Concerns 

In the previous discussion, we have assumed accretion given by the Bondi/Hoyle rate 
with radiation feedback. Each seed black hole will accrete more rapidly with time, and 
there can be an epoch when the seed black holes reach the Eddington limit (for the relevant 
opacity). In this case, the accretion would proceed in a different manner. 

We can solve for the time at which the Eddington limit for electron scattering is reached 
during the Bondi accretion phase by setting the accretion luminosity for a single black hole 
from Equation |31] equal to the Eddington luminosity: 

L{MBH{t)) = vFMIhC^ = Lem = KMbh (49) 

where KthI^]^ erg s"^ g^^. We thus have for the black hole mass at the Eddington limit, 

MBH,Edd, 

Mbh,em = ^ ~ 250 Me n-^Hf^ (50) 

for fiducial values of temperature parameters in Equation [2TJ The time when the first seed 
reaches the Eddington limit is obtained by equating MBH,Edd with MBH,max from Equation 
|32| neglecting the delay time, or, 

i\/r ^ i\/r Mbh,o /^cin 

vFc^ 1 - FMsHfiitEdd - to) 

or, 

K - tjFc^Mbh 
tEdd -to = ■ too, (52) 

A _^ IVlBHfi 

where in the final step we have neglected the second term in the numerator that is numerically 
small compared to K. Thus, formally we would have for Mbh < ^ ' 



BH = J^ ^^^BH 



MnH = FM^ 
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M = ^^^.0 

'''' l-F{t- tacc)MBHfi 

L{t) = vFc^MIj,, (53) 



and for Mbh > ^: 



M, 



KM 



BH 



BH 



rjc^ 



Mbh = —^r^e^ 

L = KMbh. (54) 

In practice, tEdd ~ ^oo is so long even under conditions of high density, Uism ~ 10^ cm~^, 
that it is unhkely that any of the seed black holes would reach the Eddington limit before 
the ambient conditions had changed to render the Eddington limit even further out of reach. 
For the conditions we envisage, the Eddington limit has no practical influence on the growth 
and radiation from the ensemble of stellar mass black holes. 



5. Discussion and Conclusions 

A plethora of stellar-mass black holes accreting in young galaxies could be significant 
source of radiant energy. Unlike the feedback from stars and supernovae that are one- 
time, rather short term events in the history of a galaxy, individual stellar-mass black holes 
accumulate and continue to radiate as long as the conditions promote accretion. A principal 
factor is the ambient density of the ISM that must remain sufficiently high. We have shown 
that for predicted ambient densities in the range 10^ — 10^ cm~'^, the condition for appreciable 
accretion is met in the context of Bondi/Hoyle accretion limited by radiative feedback. In 
this case, there is potentially a significant amount of accretion-radiated energy available. 

We derive the total energy production of stellar mass black holes accreting from a 
high-density primordial ISM using the classical Bondi spherically symmetri c accretion rate 



scaled with the value for the mean accretion rate affected by feedback given in lPark fc Ricotti 



(J201ll ). Since stellar mass black holes will be distributed in proportion to star formation, their 
energy production will be emitted more uniformly over the galaxy than a single supermassive 
black hole. This may have implications for early galaxy evolution as large scale galaxies start 
to form during the era of reionization. We also show that given the estimated values used 
in our calculations, the stellar-mass black holes are not likely to reach the Eddington limit 
during the time span we have used here, a few Gyr. 
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For our analytic calculations, the ISM density, riism, has been taken to be constant in 
space and over the time spans of interest. In reality, the density will vary througout the 
galaxy and as the galaxy evolves over time, the density will decrease. The motion of the 
black holes with respect to the background gas could also be an important limiting factor 
as the motion approaches or exceeds supersonic. This motion could affect the nature of 
the feedback process itself, but, to the best of our knowledge, this has not been explored 
in any detail. Due to these effects, our calculations may be more optimistic than a more 
realistic case that accounts for density fluctuations and evolution and for black hole motion. 
Some account of these factors might be taken by replacing the density of the ISM in our 
formulation, riig^n with ffiiism where fj is a filling factor corresponding to a given density. 
The ambient density may be affected by whether star formation is dominated by mergers 
or by cold inflow at the epochs of interest (jSomerville et al.ll2008l ). We have, however, given 
a general framework by which the accumulating number of stellar-mass black holes, their 
luminosity at a given epoch, and their accumulated feedback energy could be incorporated in 
a numerical simulation that followed the fluctuations and density evolution more realistically. 
A simulation that included the distributed feedback effects of accumulating, accreting stellar- 
mass black holes might give significantly different results than one that only included stars, 
supernovae and one or a few intermediate-mass or supermassive black holes. 



As we were finishing work on this paper, the paper by iMirabel et al.l (l201l[ ) appeared 
that considers the possible role of stellar-mass black holes in high-mass X-ray binaries in 
young galaxies. This scenario has the advantage that the accretion is driven by mass transfer 
from a companion star and does not depend on the vagaries of the ISM. On the other hand, 
this proposition does depe nd on the binary frac tion of massive stars in the early Universe, a 
rather uncertain quantity. iMirabel et al.l (120111 ) seem to assume that all stellar-mass black 
holes formed in these early epochs are in binary syst ems. We note that these high- mass 
X-ray binaries live for a short time (perhaps 0.02 Gyr, IMirabel et al.l 1201 if ) and hence each 
binary system is a delta-function contribution on the time scales we consider here. When 
the secondary dies, the primary black hole (or a secondary one if it forms) would then be 
subject to the sort of radiative-feedback limited Bondi accretion we consider here. Both 
single and binary black holes should be considered in complete models of the evolution of 
young galaxies. 

While we have concentrated on epochs ~ Gyr, the existence and impact of stellar-mass 
black hole accretion and feedback may be pertinent to the earliest phases of star and galaxy 
formation if, for instance, turbulence allows the formation of a broader initial mass function 
and hence lower mass stars already in the first dark matter halos where star formation is 
thought to have first begun (I Clark et al.l 1201 ll : iPrieto et al.l 1201 if ) . Whereas a single very 
massive star explosion by pair-instability releasing ~ 10^^ ergs is likely to unbind a mini-halo. 
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the explosion of less massive "normal" supernovae and the accretion and feedback of stellar- 
mass black holes may yield a different behavior than many current simulations envisage. 
Supernovae may blow fountains rather than totally disrupting the baryonic content of the 
halos and feedback from accreting black holes may alter state of ionization of the gas. Both 
types of events may drive further turbulence. If stellar-mass black holes are present from the 
first epochs of star formation, then their role must be considered in the subsequent mergers 
that build more massive galaxies. Aside from the effects of accretion and feedback, any 
stellar-mass black holes that form and remain in mini halos may themselves merge through 
dynamical friction and help to promote the growth of larger-mass black holes during the era 
of rapid galaxy merging. The tendency to move supersonically with respect to the ambient 
gas that will tend to limit the accretion may promote this dynamical friction and merging. 

An important aspect that we have not yet explored in depth is the observational con- 
sequence of our model. This might be addressed by a simulation incorporating stellar-mass 
black hole feedback. A related issue might be to discriminate the signal of a host of stellar- 
mass black holes from single, larger-mass black holes. One means to do this might be to make 
both radio and X-ray observations of young galaxies. There appears to be a nearly universal 
relation between the radio and X-ray luminosity of black holes such that L^adio ot L^'^ , where 
both luminosities are presumed to scale monotonically with black hole mass. If this holds 
for the conditions we address here, then an ensemble of stellar-mass black holes should have 
a larger radio fl ux for a given X -ray flux than a single black hole of the same total mass. 



In this context, iJia et al.l (|2011[ ) have examined the evidence for the growth of black holes 
in local starburst galaxies. They flnd an anomalous increase in X-ray flux that they argue 
could be associated with the growth of a black hole of lO^M© by accretion. We suggest there 
might be an ensemble of 10^ smaller-mass black holes that might provide the same X-ray 
flus, in which case the radio luminosity would be relatively large. 

Combining Eqns ^U\ and [21] we estimate the accretion rate for flducial parameters to 
be ~ 10~^ Mq yr~^. This accretion ra te corresponds roughly to that need ed to induce 
the accretion disk limit-cycle instability ( ICannizzo et al.lll982l : iLin et al.lll985f) and thus to 



prod uce a black hole X-ray nova analogous to AO620-00 and related events ( iTanaka &: Lewin 



19951 ). The X-ray flux from these systems might thus come in flares near the Eddington 
limit lasting for months with quiescent periods of years to decades. In an active star-forming 
galaxy at a red shift of about 2 with ~ 10^ stellar mass black holes, ~ 10^ of them might 
be in outburst at any given time. Any such bursts may be redshifted into bands that are 
heavily extincted and hence difficult to observe directly, but this aspect is worth considering 
more carefully. 



Kormendy et al.l (120111 ) have argued that black holes do not correlate with disks and 
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that they correlate little, if at all, with pseudobulges. It would be interesting to consider 
whether or not stellar-mass black holes affect the early evolution of disk-grown pseudobulges. 

There is currently a problem understanding the rapid early growth of dust in young 
galaxies. Perhaps the ensemble of stellar-mass black holes contributes to dust formation 
by, for instance, providing numerous and prolonged concentrations of ISM density at the 
boundary of radiative feedback regions. 

Finally we note that our estimates here suggest that a contemporary stellar mass black 
hole that found itself in a dense clump of molecular material, n ~ 10^ cm~^, might be 
observably luminous ~ 10^^ erg s~^ for fiducial parameters at this density. Dense clumps 
like this are estimated to fill a volume in the Galaxy of about 4 x 10^ pc^ (N. J. Evans, 
private communication). With a volume within 100 pc of the Galactic plane of ~ 10^^ pc^, 
this represents a filling factor for dense molecular clumps of ~ 4 x 10~^. With estimates of 
10^ — 10^ black holes in the Galaxy, the estimated number of black holes within such dense 
clumps ranges from negligible to a few. 

We thank Milos Milosavljevic, Pawan Kumar, Volker Bromm, Jenny Greene, and Neal 
Evans for useful discussions. Special thanks go to Tracey Bennett who helped with the 
literature search. VJ is grateful for an award from the College of Natural Sciences and for 
support from the Undergraduate Excellence fund of the Department of Astronomy. This 
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Fig. 1. — Schematic diagram of the growth of a black hole from a seed by feedback-limited 
Bondi/Hoyle accretion in an early galaxy. 
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Fig. 2. — The mass of a stellar-mass black hole as a function of the time the black hole 
accretes for an initial mass of 10 M© in an early galaxy with star formation rate 10 M© yr~^ 
for various values of the ambient density, riism- Other parameters are set to fiducial values 

(see text). 
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Fig. 3. — The number of stellar-mass black holes as a function of time in an early galaxy 
with star formation rate 10 M© yr~^ for various values of the initial black hole mass, Mbh,o- 
Other parameters are set to fiducial values (see text). 
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Fig. 4. — The luminosity of stellar-mass black holes as a function of time in an early galaxy 
with star formation rate 10 Mq yr~^ for various values of the ambient density, nism- Other 
parameters are set to fiducial values (see text). 
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Fig. 5. — The luminosity of stellar-mass black holes as a function of time in an early galaxy 
with star formation rate 10 M0 yr~^ for various values of the velocity parameter f^ = 
{v^ + cj.)/cl and riism = 10^ cm~^. Other parameters are set to fiducial values (see text). 
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Fig. 6. — The total integrated energy of stellar-mass black holes as a function of time in an 
early galaxy with star formation rate 10 M© yr~^ for various values of the ambient density, 
nism- Other parameters are set to fiducial values (see text). The curves are truncated for t 
approaching too from below due to graphics resolution. 
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Fig. 7. — The total integrated energy of stellar-mass black holes as a function of time in an 
early galaxy with star formation rate 10 M© yr~^ for various values of the velocity parameter 
fv = (v^-|-c^)/Cg and riism = 10^ cm~^. Other parameters are set to fiducial values (see text). 
For larger values of f^, the parameter too, the time for an initial black hole to grow to infinity, 
is large and the total energy has an asymptotic behavior (see text). 
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Fig. 8. — The total integrated energy of stellar- mass black holes as a function of the pa- 
rameter too, the time for an initial black hole to grow to infinity, for various fixed epochs, t, 
in an early galaxy with star formation rate 10 Mq yr~^. The curves are truncated for too 
approaching t from above due to graphics resolution. 



